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Abstract 

The work of Davis [I], which imports the concept of total-variation- 
diminution (TVD) into non-upwinded, Lax-Wendroff type schemes, is reformulated 

in a way which is easier to analyze. The analysis reveals a class of TVD 
schemes not observed by Davis. 
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Introduction 


Davis [1] has recently sought to show that concepts developed in the 
context of upwind-differencing schemes, such as flux limitation and total 
variation diminution, can be applied to more traditional, Lax-Wendrof f-like 
algorithms, to provide a rational artificial viscosity. His technique was to 
take a particular class of upwind TVD schemes in a form analyzed by Sweby [2], 
and then to modify them in such a way as to make them independent of the wave 

direction. 

Our aim in this note is to show that results similar to, but somewhat 
stronger than, those of Davis can be obtained more simply and directly. Only 
the case of one-dimensional linear advection will be treated. An extension to 
two-dimensional nonlinear systems is given in [1]. 


Analysis 

To solve the equation 

+ ai^ =0 ^ ^ 

on a regular rectangular mesh x = iAx, t = nAt with mesh proportions given 
by the Courant number v = aAt/Ax, consider schemes of the form below, in 
which Au 1+1/2 = u" + 1 - u" 
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n+l n 


u t - u 4 = - V 2 v(l + v)Au ± _ I/2 v(l - v)Au i+ \j 


- V 2 M(1 - M)(l - Q._ i / 2 >Au._ 1/2 


+ V 2 |v|d - |v|)(l - Q i+ 1/2 )Au i+ 1/2 . 


( 2 ) 


Here the first line represents the usual Lax-Wendroff scheme, and the 
other terms represent a conservative dissipation, in which + y depends on 
three consecutive gradients, Au._ ^ , Au ±+ y , Au i+3/2 . In fact, since 

Qi+ l/ 2 is dimensionless, it can only depend on the ratios of those gradients, 
and so we write it as 


Au , 1 , 

Q - Q( 1 / ; , 

1+ V 2 4 V l / 2 


4u 1+3/2 j 
A “i+ 1/2 


(3) 


or, more concisely 


Q 


i+ V 2 



r i + v 2 


). 
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The particular form chosen for the scheme will justify itself when it 

turns out to be very simple to analyze. The factors | v | ( 1 - |v|) 

multiplying the dissipative terms can be motivated by noting that the basic 
Lax-Wendroff scheme is exact and needs no modification if |v| = 0 or 1. We 
shall show that Q can be chosen in such a way as to ensure that (2) 
total-variation diminishing, in the sense of Harten [3]. For this purpose, 

Davis made use of the Harten-Sweby lemma which requires the scheme to assume 
the form 

u+1 n _ 

U i ■ "l - C i- l/ 2 4 V l / 2 + »!+ V 2 4u l+ 1/2 • (5 > 






Jta gu a. 




which can be rewritten 


U i = Cu i-1 + (1 " C ” D)u i + Du i+1* (6) 

Evidently , if the weights C, D, 1 - C - D in (6) are all positive, u™ +1 

will be bounded by the greatest and least of u? , , u?, u? and this is 

l-l i l+l 

sufficient to ensure that the total variation of u n+ * is less than that of 
u n . A somewhat stronger constraint, which we shall employ here, is to 
require, when v> > 0, that D = 0 and 0 < C < 1 (and when v < 0, that 
C = 0 and 0 < D < 1). This condition specifies that must be bounded 

by the data in the "upwind^ interval. It is rather surprising that this 
upwind constraint can be mot by a non-upwind scheme. 

To show that it can, consider first the case v > 0, and rewrite (1) as 


n+1 n 

u. - u. = 
i i 


- vAu. i 

i- J 


■ A 


1 - v)AUj_ jy - 1/2 \>(1 - v)Au i+ 1 


h 


l /2V(l - v)(l - Qj_ i/ 2 )Au ± _ y 2 + V2 v(l - v)(l - Q i+ \/ 2 )Au i+ i /2 


= - vA 


u i _ l / 2 + V 2 V( 1 - v)Q._ 1/2 Au._ 1/2 - 1/2 V( 1 - v)Q. + 1/2 Au. + i /2 


v[l - V 2 (1 - V)Q._ 1 / + V 2 (1 - v)Q i+ j, /r' + j ]Au._ y 


(7) 


which is of the form ( S'i with D = 0, and 

C - v[l - V 2 O - V)Q X _ 1/2 + 1/2 0 - v)Q 1+ 1/2 /r- + 1/2 J. 


( 8 ) 
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The condition that C is positive yields 


^i- V 2 ^i+ V 2 ^ r i + V 2 1 - v 


( 9 ) 


and the condition that C is less than one is 


^i+ V2 /,r i+ V2 V2 ^ v 


( 10 ) 


The case v < 0 follows a similar pattern, which is most clearly revealed 
by writing |v| = - v, so that 


u i +l " u i = M Au i+ i / 2 “ l v l O ~ M > Au i+ i / 2 + x !l i v l (l “ l v l) &u i-i / 2 


-V2M0 - M)(l - Qj. y 2 )^ i _ l / 2 
+ V2M0 ~ M)(i ” Q i+ 1/ 2 ) Au i+ 1/ 2 
= MAu i+ 1 / 2 - V 2 M (1 - M)Q i+ i / 2 Au i+ i / 2 

+ v 2 i V I o - M)Q i _i /2 Au ± _ i /2 

= M(i - V 2 ( 1 - M)Q i+ i/ 2 + V2O - H>W 2 /r i- 1/2 ^ Au i+V2 (n) 


which is of the form (5) with C = 0 and 


D= | v | [ 1 - V2 ( 1 - |v|)Q i+ i/ +V2O- l v l >Q i _ \j 2^1- l / 2 J* 


■>» 
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The similarity between (12) and (8) reveals that to ensure 
conditions are 

^i+ V 2 V 2 ^ r i- V 2 ^ 1 - I v | 

Q i- V 2 /r i- V 2 " Q i+ 1 / 2 < M * 

The task of devising a function Q which meets the conditions (9), (10) 
(13), (14) is greatly simplified if it is assumed that both Q and Q/r are 
always positive. In that case we have 


0 < D < 1 the 

(13) 

(14) 


Q i+ 1 / 2 < 1 - 2 |v| 

Q i+ 1/2 /r i+ l/ 2 < flr 

Q i+ l/ 2 /r i+ l/ 2 < M * 


These inequalities are precisely those which appear in the theory of flux- 
limited upwind schemes using "B-f unctions" [4] which depend only on r~ when 
v > 0, and only on r + when v < 0. B(r) is bounded as in (15), and 
B(r)/r as in (16) or (17). Here Q must bear a bounded ratio to each of its 
arguments. 

To establish a connection between the present analysis and that of Davis 
[1], equation (1) of this note should be compared with equations (3-16) and 
(3-18) of [1]. Using the present notation of r”, r + , which is slightly 
different from the notation used by Davis, it will be found that the 
algorithm, are identical if 
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Q K+ 1/2 ’ r i+ 1/2 ) = *^ r i+ 1/2 ' + ^ r i+ 1/2 ^ “ 1 ’ (18) 

where <J>( r ) is Davis' limiting function. This shows that Davis has in effect 
considered the special case of Q— functions which are 'separable' in the sense 
that they are the sum of two functions each depending on one of the variables. 
For such functions Davis establishes the TVD property (not the stronger 
property proved here) provided 

0 < *(r) < 1 (19) 

0 < ij>(r)/r < 2. (20) 

These "separable Q-functions" do not necessarily obey conditions (15) - (17), 
(Q may go down to -1.0, and Q/r may not be bounded). Thus the Q-functions 
studied here are more general, in the sense of having no special functional 
form, but more restricted, in the sense of meeting a stronger condition. 

To show examples of each type, define "minmod" to be the function which 
returns the smallest number from a list of positive arguments, but equals zero 
if any argument is negative. Then a common, though not particularly good 
limiter for upwind schemes is 4-(r) = minmod (l,r). Based on this, we could 
define the separable Q-function 

Q(r ,r + ) = minmod (l,r ) + minmod (l,r + ) -1, (21) 


or the non-separable Q-function 


Q( r >r ) = minmod (l,r ,r + ). 


( 22 : 


both of which are sketched below. 

We next turn to the question, whether any special choices of the Q- 

function .ill produce algorithms with distinguished properties, such as third- 
order accuracy. It is easily shown [4] 



Equation (22) 


I 

r 



0 

r 

1 

r + - 1 

— 4 

r + r - 1 

r 

0 



-1 

r~ - 1 

0 


Equation (21) 


that upwind schemes using a limiting function such that when 
unity 


r is close to 


♦<r) r + 

3 3 


(23) 


are third-order accurate 
seek a linear Q-function, 
the form 


in smooth regions of the flow. Analogously, we may 
to be used if r , r are both close to unity, of 


Q(r ,r + ) 


ar + b + cr + . 


1 ">/. \ 
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Substituting this expression into equation (6) produces the algorithm 


u n+1 - u? = - vAu. i. + •§■ v(l - v)Au . . + v(l - v)Au, 

li i- 72 2 1-3/2 2 i 


-v 2 


+ |v(l - v)Au 1+1/2 -| vd - v)Au._ i/2 


-|v(l - v)Au. + 1/2 -| v(l - v)Au. +3/2 


(25) 


Conditions for the accuracy of algorithms of this form were given by Roe 
[5], the general case being 


n+1 n v * 

u . - u . = -v ) y Au ( , |, . 

i i L k i+V2“ k 


(26) 


with the conditions for first, second, and third-order accuracy being 
respectively 

l Y k = 1 (27a) 

k k 


}' ky, = V 2 ( 1 + v) 


(27b) 


l k 2 Y k = 1/6(1 + v)(l + 2v) , 


(27c) 


In equation (20) the coefficients are 


Y_ 1 = V 2 (1 " v )c 


(28a) 


Y 0 = V 2 ( 1 ~ v)(b - c) 


(28b) 
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Yj = V 2 (1 - v)(a - b) + 1 


(28c) 


Y 2 = - V 2 (1 ~ v)a. 


(28d) 


Obviously (27a) is satisfied for 


any a, b, c, and (27b) is satisfied provided 


a + b + c = 1 , 


(29) 


After simplification, (27c) reduces to 


3a + b - c = ( l/3)(2v - 5). 


(30) 


Clearly, there is no way to soecifv => , . 

y co specity a, b, c independently of v so that 

(30) is satisfied, and our ambition to create a a 

create a third-order, non-upwind, TVD 

scheme is frustrated. 

in the absence of a third-order scheme, we may seek special second-order 
schemes. One possibility is to preser.e the property of the basic Lax- 

Wendroff scheme that it connects exactly any quadratic function of x. For 
any such function three consecutive differences 4„ are in arithmetic 
progression. Thus the arguments of Q are of the form 


r - 1 - e. 


r = 1 + e 


(31) 


and we should seek Q-fu„ctio„s which equai unity for these arguments. The 
dissipation terms in equation (,) *11 then vanish. For example, any linear 
function such as (24) in which a - c and a ♦ b ♦ c - , win have thls 
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property. It is worth remarking that Q = ty 2 (r + r + ) generates a fourth- 
order dissipation. Of course, linear Q-functions cannot be used for large 
values of the arguments because they exceed the bounds placed on them by 
conditions (15) - (17). However, for arguments close to unity, linear Q- 
f unctions can be used, although they must be replaced for other arguments. As 
an example, consider the readily-evaluated function 


Q = minmod (2, 2r~ , 2r + , V 2 (r~ + r + )) 


(32) 


whose behavior is displayed in the sketch. 


+ 



This function will convect exactly any quadratic data, provided the gradient 
ratios, r , r do not lie outside the triangular region. That triangular 
region can be made much larger if advantage is taken of the way the limits in 
Theorem I depend on v. For example, if v = V 2 an acceptable Q-function is 


Q = minmod(4, 4r , 4r + , V 2 (r + r + )) 


( 33 ) 


and the v- rtices of triangle move to (1,7) and (7,1), 


-li- 


lt is this freedom to match the dissipation to a unique Courant number 
that wili be lost when dealing with systers of equations. Davis { 1 J suggests 
matching the scheme to the Courant number of the fastest wave. One might also 
consider matching the scheme to the strongest wave, by some such device as the 
following. Let I 3 be the wavespeeds, and ctj , the 
amplitudes, of the waves occurring in the interval i + I /2 • The values of 
these quantities for the Euler equations are well known; 


0 a, X~ — u, X — u + a 


(34) 


a — [Ap - paAuj, a = — [a 2 Ap - Apj, a = — - [a p +paAuj. (35) 


1 2 a 2 


2 2 
a 


3 2 a 2 


Define the mean wave speed X to be 


x * _ a l X 1 + a 2 X 2 + a 3 X 3 
2 2 ~2 
a l + a 2 + a 3 


(36) 


Since X is a convex combination of u - a, u, u + a, it must lie between 
u - a and u + a. It is actually given by the formula 


X* = U + ^ 2p-i. 2 A - P Au / q 7 \ 

o 99 9 ^ * \ v3 / J 

(Ap) 4- p a (Au) + 2(Ap - a Z Ap) 

In the case of two acoustic waves having equal (or opposite) amplitudes 
equation (31) returns the particle speed, which is not very useful. Perhaps 
more useful, and more in the spirit of the enterprise (since the dissipative 
terms only require the absolute value of v) , would be to compute an r.m.s. 
wavespeed, given by 
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x 2 


2.2 22 2 ,2 
a l X 1 + a 2 X 2 + °3 S 
2 2 2 
a l + a 2 + a 3 


2 2 
= u + a 


(Ap) + 4puApAu 


2 2 /a >.2 
p a (Au) 


(Ap) ^ + p 2 a 2 (Au)^ + 2 [Ap - 


(38) 


Ap 


The term which appears in the denominator of the fractions in (37), (38) 
measures the total strengch of the disturbance in the cell. It is therefore a 
candidate for the quantity whose ratios in consecutive cells will serve to 
define r~, r + (see Davis [11, equation (4.13)). All these possibilities 
however, require extensive numerical testing. 


Conclusions 

The inspiration of Davis [1] to introduce the TVD concept into non-upwind 
algorithms has been reformulated in a way which permits more general results 
to be deduced by rather simpler arguments. A new class of non-separable 
limiters emerge from the analysis. Preliminary numerical experiments have not 
shown any striking advantage to these limiters, but the simplified analysis 
should be of advantage when attempting to extend these ideas to provide 
viscosities for symmetrical algorithms other than Lax-Wendrof f . 
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